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ABSTRACT 


In the present study, boiling of highly wetting dielectric fluid has been investigated 
in an oscillating fluid environment. A piston is designed to create oscillations in the fluid 
and over the heated platinum wire. Because of their low surface tension, these liquids 
require very high superheat to initiate nucleate boiling. It is expected that the amount of 
necessary temperature overshoot for the onset of nucleate boiling, can be decreased with 
oscillation in the fluid. The oscillation can remove the bubbles, which are forming in the 
nucleation sites as soon as they start growing on the outer surface. This increases 
efficiency of nucleation sites, which are very scarce. 

All of the oscillation amplitudes and frequencies, tested here, changed the boiling 
curve of highly wetting dielectric fluid, so that the apparent temperature overshoot has 
decreased. Remarkably at some oscillation amplitude and frequencies the superheat is 
almost vanished. The effects of the amplitudes and frequencies on the boiling curve varied 
because of the present bubble size and growth rate, which depend on the size and shape of 
the nucleation sites. 
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I. INTRODUCTION 


A. BACKGROUND 

The developments in electronics and computer technology have led to higher 
packaging densities of computer chips. As a result of these developments the power 
dissipation of the computer chips has increased as well. In order to achieve lower failure 
rates during operation, longer life times and lower maintenance costs, electronic 
equipment must be cooled. If cooling is not sufficient, the reliability of electronic 
equipment will decrease. Because of inefficient cooling, component temperatures may 
exceed the limiting values, which is of the order of 85 C for most of electronic 
components. 

For quite a long time natural and forced convection cooling with air have been 
implemented as solutions to the electronics cooling problem. The disadvantages of these 
methods and increasing need for better cooling with super computers has led the way to 
immersion cooling with liquids. Immersion cooling is a method to deal with the levels of 
heat flux and volumetric heat generation, which are characteristics of high performance 
computers. This cooling technique has two distinctive advantages. First liquids as coolants 
are more effective due to higher heat capacities than gasses. Secondly cooling with liquids 
has the capability of two phase boiling heat transfer. During boiling, latent heat can extract 
enormous amounts of heat from electronic components. 

Two phase heat transfer in immersion cooling can provide higher heat fluxes and 
more uniform temperature distributions for an array of computer chips. Chu and Agonafer 
[Ref 1] report that water can not be used in spite of its suitable heat transfer properties in 
direct immersion cooling. The best choice is dielectric fluids like fluorocarbons, which are 
also compatible with semiconductors. 

The research on electronic cooling with dielectric liquids like fluorinerts (FC series 
of 3M corporation) is still going on. The FC-72 used in this work series is non-corrosive, 
dielectric and has a relatively low saturation temperature of 56 C. Liquid evaporation 



cooling can handle heat fluxes between 5 - 120 W/cm^ [Ref. 2], The peak heat flux for 
fluorinerts of 3M corporation [Ref 3] for pool boiling is about 12 - 15 W/cm^. 

Boiling with fluorinerts like FC-72 has some problems. FC-72 has very low surface 
tension, which means almost zero wetting angle. This fluid can fill in even the smallest pits 
and cavities, eliminating most of the nucleation sites on the surface, which are necessary 
for the onset of boiling. As a result of the absence of nucleation sites, the surface 
temperature departs from the saturation temperature significantly during boiling, which 
causes wall superheat. The surface temperature excursion during boiling due to the heat 
flux is the core of the problem. Absence of nucleation sites delays nucleate boiling and 
causes higher surface temperatures. 

Two techniques are under investigation to reduce the amount of excessive 
superheat required to initiate nucleation boiling for practicable applications of fluorinated 
hydrocarbons as dielectric cooling liquids in two phase immersion cooling. First surface 
properties and surface treatment can be selected such that the number of available 
nucleation sites is relatively higher. Second increasing the physical activities on the hot 
surface is expected to reduce the surface tension effects and enhance nucleate boiling. 

Surface materials and treatments, which provide more nucleation sites in the right 
range, reduce wall superheat and increase critical heat flux. Unfortunately with today’s 
technology it is not very likely to get more nucleation sites in the right size range without 
damaging the computer chip itself by surface treatment. 

Instead of improving surface properties, it is possible to decrease the amount of 
superheat, by enhancing nucleate boiling by means of physical modifications. Different 
methods of bubble pumping to decrease the effects of surface tension and detach the 
bubbles forming on the hot surface have been investigated. 

In the present study, a platinum wire is calibrated to read temperatures in pool 
boiling of FC-72. A sinusoidal motion is created by a piston in the fluid. It is expected that 
the motion of the fluid will detach the bubbles, which form on the platinum wire surface 
and decrease the amount of superheat necessary to initiate boiling. 



B. PREVIOUS WORK 

You [Ref. 4] has conducted tests which have shown that the diameter of the 
bubbles required for the onset of nucleate boiling on highly wetting dielectric fluids like 
FC-72, is 0.05-0. Igm. 

You, Bar-Cohen and Simon [Ref 5] conducted experiments with cylindrical heater 
surfaces immersed in FC-72 and R-113, saturated at one atmosphere pressure. They have 
observed the effects of fluid and surface properties and materials on boiling regimes and 
the onset of nucleate boiling. They have found that incipience wall superheat reduced 
about 25 C and nucleate boiling superheat decreased about 6 C in FC-72 compared to R- 
113. They have concluded that the fluid properties are more effective on boiling curve 
than surface material and FC-72 is more advantageous for cooling than R-113. 

You, Simon and Bar-Cohen [Ref 6], studied enhanced surfaces with more 
available nucleation sites in immersion cooling with boiling. The study showed that a 
particle layering technique, which could be applied to computer chips without damage or 
stress, lowers the wall superheat in boiling and increases critical heat flux compared to the 
untreated surface, 

Marcellus [Ref 7], conducted experiments to show that forced flow can increase 
the critical heat flux value by rapid removal of vapor bubbles from the surface. A peak 
heat flux value of 78 W/cm^ with FC-75 has been achieved, when the fluid velocity was 7 
ft/sec. 

Recently, Kelleher, Egger Joshi and Lloyd [Ref 8] conducted experiments on a 
vertical array of five horizontal platinum wires. They have observed that a temperature 
overshoot is necessary to initiate nucleate boiling in highly wetting dielectric liquids. They 
have conducted some experiments to determine the effect of the bubbles coming from the 
lower wires on the incipience of nucleation boiling on the surface of the upper wires. Their 
results show that the bubbles from the upstream source reduced the temperature 
overshoot on a downstream surface. They have used FC-72 as the dielectric fluid in their 
experiments. 





C. OBJECTIVES 

The purpose of this experimental study is to achieve superheat reduction by 
forcing bubbles to depart from the hot surface with the help of fluid motion which 
overcomes surface tension effects. Departing bubbles provide the nucleation sites for new 
bubble formation and give way to nucleate boiling. The fluid motion is created in FC-72 by 
piston displacement. Different piston stroke lengths and frequencies should be tested to 
see the temperature excursion on the wall. 




n. EXPERIMENTAL APPARATUS 


A. DESCRIPTION OF COMPONENTS (Figure 1) 

The experimental setup consists of four main parts: the main chamber; the front 
chamber; the adjustable stroke piston, and DC motor. The main chamber has been used in 
the previous experiments to investigate the boiling of FC-72 under different circumstances 
and is modified to be used for this study. The goals of this experiment required some 
modifications, which are the design of the adjustable stroke piston, DC motor and the 
front chamber. The front chamber regulates the motion in the fluid, which is caused by the 
piston, and transmits the motion to the main chamber. Figures 2a and 2b show a different 
view of experimental apparatus and its schematic drawing. 

The main and the front chambers hold the dielectric fluid FC-72. To prevent lose 
of FC-72 during degassing, sealing is very important. At every connection point, RTV 
Silicone Adhesive Type 1 is used as a sealant. At the clearance between the piston seal and 
the front chamber, teflon tape is utilized to prevent leakage. 

1. Main Chamber (Figure 3) 

The main chamber is a Plexiglas box with an aluminum cover on the top. The inner 
dimensions of the main chamber are 5.9964 x 2.3776 x 5.9843 in. (Figure 4a). This 
measurement is necessary to calculate the displacement of the fluid in the chamber and the 
amplitude of the motion during the experiments with the piston running at a given stroke 
length. 

The board, on which the 0.05 mm platinum wire is attached with an adjustable 
spring system, can be placed in and removed from the chamber by means of two slots on 
the side walls. In order to connect the heaters and the platinum wire to their power 
supplies and the thermocouples to the data acquisition system, a tube (Figure 3) at the 
back wall of the chamber is used. All the wiring necessary for the equipment in the main 
chamber exits via this tube and goes either to the connection board or to the data 
acquisition system 



In order to reduce buoyancy induced convection from the bulk heaters one strip 
heater was used to maintain the FC-72 close to saturation temperature. When the 
experiment were conducted close to the saturation temperature, heat losses were critical. 
Because of this the main and the front chambers both are covered by two layers of 
insulation to decrease the heat losses. 

Four thermocouples in the main chamber read the bulk temperature of the FC-72. 
The passage between the front chamber and the main chamber is an oval opening and the 
front chamber side of the opening is covered with a flow regulator, to provide uniform 
flow from the front chamber to the main chamber. 

The main chamber contains three important parts: the condenser; the bulk heaters; 
and the board with the adjustable spring system. 

a. Condenser (Aluminum Cover Plate) (Figures 4a, 4b and 5) 

The main chamber is covered by an aluminum plate, which condenses FC- 
72 vapor and returns it to the chamber. The aluminum plate is fixed to the chamber by 8 
screws, which apply required pressure on the 0-ring at the same time. As it has been 
mentioned, for additional insulation this connection part is covered with RTV sealant. The 
lower surface of the cover plate has 4 thermocouples embedded in a cross shaped groove, 
(Figure 4a) which read the condenser surface temperature. 

On the upper surface of the condenser (Figures 4b and 5), there are three 
thermo-electric coolers connected in parallel, which are manufactured by MELCOR 
THERMOELECTRICS [Ref 9]. They keep the condenser temperature below saturation. 
These coolers are adjustable to certain range of temperature differences. They are 
controlled by current, because the voltage drop across the thermo-electric coolers changes 
non-linearly with temperature. The thermo-electric cooler in the middle is connected to a 
toggle switch, which is another alternative to control the condenser temperature.(Figure 6) 
CP 1.4-127-045L type thermo-electric cooler have been used. These have 
an operating point of; I = 8.5 Amps, V = 15.4 Volts, and the dimensions are 40 x 40 x 3.3 
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mm. The CP series of this type of cooler provides larger heat pumping applications for 
higher currents. The thermo-electric coolers are operated at 5 A in parallel. 

To keep the cold side surface at about 20 C, the hot surface is equipped 
with an aluminum finned heat sink and air is blown between the fins by a fan. This helps to 
provide the lower and the upper limits of the temperature difference on the cold and hot 
sides as low as possible. 

b. Heaters 

Three strip heaters are used to provide the required bulk heating during 
degassing and experimentation. The heaters are fixed to a very thin aluminum plate which 
is at the bottom of the main chamber. There are two main purposes of heating in this 
experiment. Firstly, the air dissolves up to 48% [Ref 10] by volume in FC-72, that is why 
before each experiment the fluid must be degassed by boiling. Secondly, during the 
experiment the fluid must be kept close the saturation temperature, which is 56 C for FC- 
72. 

The strip heaters used in this experiment are manufactured by MINCO 
company. (HK5224R28 .9L 36A 9049). The heater dimensions are 125x11 mm. The 
heaters are in parallel connection and all three heaters are used at 2.2 A and 20 V during 
degassing and just one heater is used during the experiment at 0.6 A and 20 V. For that 
reason two of the heaters are connected to the power supply with toggle switches. (Figure 
6) The heater stays on during the experiment is the one, which causes the least disturbance 
in the platinum wire area. 

c. The Board and Adjustable Spring Assembly (Figure 7) 

The 0.05 mm diameter platinum wire and the spring, which keeps the wire 
in tension against thermal expansion, are mounted on opposite sides of a Plexiglas board. 
As can be seen from Figure 7, by the adjustable spring assembly, it is possible to solder the 
platinum wire to the brackets under no tension and then apply required spring 
displacement gradually without breaking the wire by the adjustable spring assembly. 



2. Front Chamber (Figure 3) 

The front chamber is added to the main chamber to absorb the irregularities of the 
fluid displacement, which are caused by the piston motion and transmitted to the flow in 
the main chamber. For this purpose a spongy and permeable substance is placed in the 
opening between the two chambers as a flow dampener. 

3. DC Motor and Piston Assembly (Figures 1 and 2) 

In this experiment, a piston assembly has been designed to create a sinusoidal 
motion in the fluid and therefore over the platinum wire, to remove the bubbles forming on 
the wire. Under these circumstances, it is expected that the new bubbles will form on the 
surface and they will be removed by the next wave of the sinusoidal motion. This should 
lead to nucleate boiling with less temperature overshoot from the saturation point. 

The frequency and the amplitude of the sinusoidal motion, which provide lower 
temperature overshoot from the saturation point, are to be determined. Because of this the 
variable piston stroke and frequency are needed, to investigate the possible frequency and 
amplitude ranges. 

The variable frequency is achieved by a 27 volt, 1 amp DC motor. The rotational 
speed is measured by a micro switch tachometer. The variable stroke is achieved by a slot 
(Figure 1) on the rotating disk connected to the DC motor and its fixing mechanism. The 
piston has a 0.3762 in. diameter and runs through a replaceable Plexiglas part connected 
to the front chamber with four screws. In case, the amount of displacement with maximum 
piston stroke does not suffice to get satisfactory results and larger displacements are 
required, it is possible to replace this Plexiglas piece with an appropriate size piston. 

B. INSTRUMENTATION (Figure 8) 

The main goal of this experiment is to observe the boiling curve of FC-72 by 
obtaining the surface temperature of the platinum wire at various surface heat flux values. 
For that reason a platinum wire is calibrated to temperature vs. resistance values and 
heated electrically. The platinum wire is connected to a precision resistor in series. By this 
setup every time the voltage drop across the precision resistor is known, the current in the 




line can be calculated by a simple division. Next the voltage drop across the platinum wire 
is measured and the resistance of the wire is calculated by dividing the voltage drop by the 
current. With the known resistance value, the surface temperature of the platinum wire can 
obtained from the calibration curve. Since the voltage drop and the current are measured, 
their product yields the heat transfer rate and division of the heat transfer rate by the 
surface area of the platinum wire gives the surface heat flux. 

The bulk temperature of the fluid and the condenser temperature are measured by 
the data acquisition system. Some additional instrumentation is required to obtain the 
frequency of the DC motor and the piston. Figure 6 is the wiring diagram of power 
supplies and data acquisition system without thermocouples. 

1. Platinum Wire 

The platinum wire is a product of Goodfellow Advance Materials. It has 99.9% 
purity and 0,05 mm diameter. The length of the platinum wire in this experiment is 4.011 
in. The Scanning Electron Microscope investigation of the sample wire from the same 
spool of the wire used in the experiment yielded 55pm diameter. The result of the 
platinum wire calibration is a linear temperature vs. resistance relationship. Moving from 
calibration bath to the experiment apparatus and especially inserting the board into the 
main chamber, requires very careful handling of such thin wire. 

The resistance of the precision resistor, which is in series with the platinum wire, is 
m ± 1%. The precision resistor are manufactured by Dale company and have 25 watt 
power rating, 

2. Power Supplies 

In this set of experiments, four DC power supplies are used for the platinum wire 
and precision resistor in series, thermo-electric coolers in parallel, heaters in parallel and 
DC motor. These are: 

HP 6289A DC Power Supply 0^0 V, 0-1.5 A with HP 59501B Isolated DAC/ 
Power Supply Programmer.( HP 6289A provided voltage to the platinum 
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wire and precision resistor and the voltage value is controlled by HP 
5950IB and the computer program.) 

HP 6286A DC Power Supply 0-20 V, 0-10 A (to power the thermo-electric 
coolers.) 

HP 6269B DC Power Supply 0-40 V, 0-50 A (to power heaters.) 

HP 6289A DC Power Supply 0-40 V, 0-1.5 A (to power DC motor.) 

3. Data Acquisition Unit and Communication Board 

All the temperature and voltage measurements are done by HP 3497A Data 
Acquisition system via Unitek desktop computer with Pentium processor. An interface 
adapter provides communications between the computer and the acquisition system. The 
interface adapter is a National Instruments, GPIB-PCII/IIA with NI-488.2 software for 
DOS and Windows (part # 776092-01). The same interface adapter also allows the 
implementation of HP 5950IB Power Supply Programmer. 

The cards, which are used for temperature and voltage readings for data 
acquisition system, are as follows; 

1. Thermocouple card 44422A 20 channel T-couple Acquisition. 

2. Voltage card 44421A 20 channel Guarded Acquisition. 

4. Frequency Measurements 

Since the frequency of DC motor is very low, DC motor period is measured by 
5001 Universal Counter-Timer and the frequency is obtained from period. 







COOLLKS 

'condense#^ f „ . 


Figure 2a. Main Parts of Experimental Setup 
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Figure 2b. 


Main Parts of Experimental Setup 











































Figure 3. Main Chamber and Front Chamber 















Figure 4a. Condenser and thermocouples 


















THICKNESS OF THE ALIMINUM PLATE = 6 nun 

THICKNESS OF THERMO-ELECTRIC COOLERS = 3.3 mm 


Figure 5. 


Schematic View of Condenser and Coolers 


































Figure 6, 


Data Acquistion Connections and Wiring Diagram of Power Supplies 
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m. EXPERIMENTAL PROCEDURE 


A. SYSTEM PREPARATIONS 

Prior to each experimental run, the following preparation steps must be completed. 

1. Check the fluid level in the main chamber and the tube. The level, in the main 
chamber, must be above the platinum wire, but should not exceed the top of Plexiglas 
board, which holds the platinum wire and the spring mechanism. The level, in the tube 
must be even with the level in the main chamber. 

2. Adjust appropriate piston stroke length before experiment. 

3, Turn on the HP 3497A Data Acquisition System. Warm up time for HP 3497A 
Data Acquisition System is approximately 1 hour.[Ref 11 pp. 648] 

4, The chambers must be insulated during degassing and experiment. 

B. DEGASSING 

Once system preparations are completed, the following steps are followed to 
accomplish degassing of FC-72, which is necessary to start an experimental run: 

1. Power all three bulk heaters using the HP 6269B DC Power Supply at 20 V and 
2.2 A. Be sure that the top two toggle switches, which are on the connection board, are 

2. Power all three thermo-electric coolers using the HP 6286A DC Power Supply 
at 4,5 V and 5 A. Be sure that the third toggle switch on the connection board, which 
controls the middle thermo-electric cooler, is on. 

3. Make it sure that the tops of thermo-electric coolers are covered by the 
aluminum finned heat sink, 

4. Plug in the fan to 110 V electric source and blow air for aluminum firmed heat 

sink. 

Watch the fluid temperature increase and the condenser temperature decrease. 
Eventually the condenser temperature must be about 20 C and fluid temperature should go 
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up to saturation. At this temperature degassing takes place with boiling. The completion 
of degassing procedure takes 1.5 to 2 hours. 

C. EXPERIMENTAL PROCEDURE 

The following steps are followed before experiment: 

1. Turn off two of the bulk heaters. 

2. Decrease the power of remaining heater to 16 V and 0.6 A. Wait till the 
temperature drops down to a thermal equilibrium point. 

3. Check piston stroke length. 

4. Turn on 5001 Universal Counter-Timer. 

5. Power DC motor by turning on HP 6289A DC Power Supply and adjust 
desirable frequency by DC motor voltage. Observe the period on the Universal Counter- 
Timer display. Stop motor after this adjustment. 

6. Turn on HP 6289A DC Power Supply to power platinum wire and the precision 
resistor. 

7. Turn on the computer and activate the interface adapter by answering Y to this 
question: DEVICEHIGH=/L: 1,51680 = C: \GPIB-PC\GPIB.COM[Y,N]? 

8. Run QBASIC, and load program EXPER3.BAS. 

9. Enter file names to store temperature and heat flux values. 

10. Run first four power settings without running the motor. 

11. At the end of fourth power setting, turn on HP 6289A DC Power Supply to 
run DC motor. 

D. DATA ACQUISITION PROGRAM 

Data acquisition is accomplished using a QBASIC program called EXPER3.BAS. 

1. The data acquisition program executes these operations: 

Sets up communications with HP 3497A and HP 6289A. 

Controls power settings incrementally by the help of HP 59501B Power 
Supply Programmer. 

Keeps time intervals between steps. 
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Reads thermocouple voltages and converts them to temperature values. 

Calculates platinum wire temperature and heat flux. 

Stores temperature and heat flux values in two separate data files. 

2. For each power setting these measurements have been taken 10 times; 

Bulk temperature (four thermocouples). 

Condenser temperature (four thermocouples). 

Ambient temperature (one thermocouple). 

Platinum wire voltage drop. 

Precision resistor voltage drop. 

3. The first five steps of power settings, the voltage increment is 0.4 V. After sixth 
step, the increment reduces down to 0.2 V. At the end of 18th reading, the program starts 
decreasing the voltage with -0.2 V till step 30 and then the decrement becomes -0.4 V at 
every step. 

4. Repeat the same steps for an other piston stroke and frequency configuration. 

E. DATA REDUCTION 

The following calculations are used to reduce the raw data obtained by 
experiments. 


1. Condenser temperature; 


Tg thru Tj are thermocouple temperatures. 


2. Fluid bulk temperature; 




(t; + t; + t; + 7)) 


T., thru T, are thermocouple temperatures. 





3. Current of platinum wire and precision resistor; 


where = 20. (Precision Resistor) 

4. Platinum wire resistance value; 

5, Platinum wire surface temperature; 


where; 


e - 50.4159 (slope of calibration curve) 

= -252.3321 (y axis intercept of calibration curve) 

6. Fluid film temperature; 




2.0 


7. Heat flux from platinum wire; 

'^Pl mre 

where = ^- 

8. Piston displacement; 

A (mm^) = — • . 

4 

9. Amplitude of oscillation in the main chamber; 

' ab 


where a and b inner length and width of the main chamber respectively. 



10. Frequency calculation, 


To evaluate the experiments, the data plotted heat flux vs. surface temperature by 
three MATLAB codes in Appendix D, Section B. Data processing has two parts; First the 
average of ten readings for each step is calculated. Second the natural convection part 
numerically calculated with recommended correlation by Kuehn and Goldstein [Ref 12] 
for natural convection around small horizontal cylinders. Then both results are plotted. 
For fluid properties, Egger’s results [Ref 13] are used. 
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rV. RESULTS AND DISCUSSION 


In this experimental study, the boiling curve of FC-72 is determined in non¬ 
oscillating bulk fluid. Then the changes in the boiling curve of FC-72 is observed when the 
bulk fluid is oscillating at different amplitudes and frequencies. All experiments are done at 
atmospheric pressure and the same platinum wire of 0.055 mm diameter is used. 

The results of these experiments are presented in the plots of heat flux vs. the 
difference of platinum wire temperature and saturation temperature in Figures 9 through 
14b. Some of the experiments are conducted for repeatability. The results of the 
experiments, which are conducted under the same conditions are in agreement within the 
uncertainty limits. 

Figure 9 is the boiling curve of FC-72 without any oscillation in the bulk fluid. The 
hysteresis between the increasing and the decreasing heat flux shows the amount of 
superheat necessary for the onset of nucleate boiling. As soon as the boiling starts the wire 
temperature decreases abruptly. Figure 10 is the boiling curve of FC-72 with the 
oscillation amplitude of 0.117 mm and oscillation frequency of 0.69 Hz. The decrease in 
the hysteresis is obvious compared to Figure 9. 

A. BOILING CURVES OF FC-72 (BOILING WITHOUT OSCILLATION) 

In order to see the effects of oscillation on the boiling regimes of FC-72, it is 
necessary to observe the boiling under non-oscillating conditions, then it is possible to 
compare the results with the oscillating cases. Figure 11 shows the results of Experiments 
4, 18 and 14, which are conducted without any oscillation. In the following discussion, the 
boiling regimes of these experiments will be compared. 

1. Natural Convection (Figures 9 and 11) 

In Experiments 4, 18 and 14. the natural convection part of the experiments occurs 
between the heat flux values of 18,700 and 86,300 W/m^ After 86,300 W/m^ is exceeded, 
nucleate boiling starts. At this heat flux value the platinum wire temperature is 99.7 C, and 
the difference between the platinum wire temperature and the saturation temperature 
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(^wire -T^,) is 43.7 C. This is the highest temperature overshoot observed during this set 
of experiments. In these experiments, the bulk temperature changes between 49 and 54.5 
C (Table la) in the natural convection phase. The bulk temperature has almost no effect 
on the natural convection curve. This observation confirms Egger’s result on the natural 
convection and the bulk temperature relation [Ref 13 pp. 23]. The data acquired in this 
part of the experiments is in agreement with the correlation suggested by Kuehn and 
Goldstein [Ref 12] for natural convection heat transfer around small diameter cylinders. 

2. Nucleate Boiling 

As soon as nucleate boiling starts, the temperature of the platinum wire drops. The 
drop in the wall temperature value changes from 11.8 to 12.6 C in different experiments 
(Table la). The platinum wire temperature just after the onset of boiling is 11.8% to 
13.5% less than the maximum superheat temperature value. 


Superheat 

Pt.Wire Temp. 

AT 

Bulk Temp. 

Heat Flux 

point 

Drop (C) 

(C) 

(C) 

(W/m") 

EXP 4 

99.6 to 87.8 

11.8 

54.5 

86,300 

EXP 18 

99.7 to 87.1 

12.6 

52.1 

86,300 

EXP 14 

99.7 to 87.3 

12.4 

49.0 

86,300 


Table la. The Comparison of Maximum Superheat Point in Three Non-Oscillating 
Experiments. 

3. Critical Heat Flux (CHF) and Entry to Film Boiling 

The location of the regime change from nucleate boiling to film boiling varies the 
most in different experiments. Different platinum wire temperatures and different critical 
heat fluxes are observed in the non-oscillating experiments for this particular point. This 
part of the boiling curve is affected most by the bulk temperature. It is observed that with 
the decreasing bulk temperature, CHF increases (Table la). In Experiment 14 (Figure 11 




and Table la), when the bulk temperature is 48.4 C, CHF is 152,000 W/m^, but in 
Experiment 18, the CHF of 132,000 W/m^ occurs at the bulk temperature of 51.8 C. At 
that time the platinum wire temperatures for Experiments 14 and 18 are 89.5 and 89.2 C 
respectively 


CHF point 

Pt.Wire Temp. (C) 

Bulk Temp. (C) 

Heat Flux (W/m^) 

EXP 4 

87.7 

54.9 

122,500 

EXP 18 

89.2 

51.8 

132,000 

EXP 14 

89.5 

48.4 

152,200 


Table lb. The Comparison of Critical Heat Flux in Three Non-Oscillating 

Experiments. CHF Increases with the Decreasing Bulk Temperature. 

B. OSCILLATION EFFECTS ON BOILING CURVE 

Three different piston strokes of 10, 15 and 20 mm, created three different 
oscillation amplitudes in the fluid surrounding the platinum wire. The amplitudes are 
0.078, 0.117 and 0.156 mm respectively. At each oscillation amplitude five frequencies are 
tested. These frequencies are from lowest to highest: 0.4, 0.56, 0.69, 0.8 and 1.5 Hz. 

1. Amplitude and Frequency Effects on the Boiling Curve 
Tables 2a and 2b represent the maximum superheat temperature and CHF values 
of the boiling curve of the non-oscillating and oscillating FC-72. Figures 12a through 14b 
are the boiling curves of the non-oscillating and oscillating FC-72. 
a. Amplitude of 0.078 mm 

The boiling curve at this amplitude and five frequencies can be seen in 
Figures 12a and 12b. The wire temperature, bulk temperature and the heat flux values can 
be seen in Table 2a for the maximum superheat point and Table 2b for CHF point. The 
drop in the maximum temperature overshoot is 8.7% for frequency of 0.4 Hz and 19.0% 
for frequency of 1.5 Hz. The heat flux at this point changes from 71,700 to 72,000 W/m^ 
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Critical heat flux decreases for three frequencies from 132,000 W/m^ to 
122,000 W/m^. For the other two frequencies CHF remains unchanged as in non¬ 
oscillating FC-72. 


Amplitude 

0.078 

0.117 

0.156 

Frequency 

(Hz) 

Wire 

(C) 

Bulk 

(C) 

Heat 

Flux 

(W/m^) 

Wire 

Temp 

(C) 

Bulk 

Temp 

(C) 

Heat 

Flux 

(W/m^) 

Wire 

(C) 

Bulk 

Temp 

(C) 

Heat 

Flux 

(W/m’) 

No Stroke 

99.7 

52.1 

86,300 

99.7 

52.1 

86,300 

99.7 

52.1 

86,300 

0.4 

95.9 

54.0 

71,700 

92.7 

53,0 

72,100 

94.7 

53.2 

72,000 

0.56 

94.7 

53.3 

71,800 

89.1 

54.8 

58,700 

93,6 

50.9 

72,000 

0.69 

94.3 

52.2 

72,000 

87.6 

54.9 

58,800 

95.3 

51.4 

79,100 

0.8 

95.7 

53.4 

71,800 

88.2 

52.3 

72,500 

95.6 

52.3 

79,100 

1.5 

91.4 

52.8 

72,200 

88.6 

53.9 

58,700 

95.1 

53.6 

71,800 


Table 2a. Superheat Values for Oscillating FC-72. For No Stroke Case the Results of 
Experiment 18 is Used. 


b. Amplitude of 0.117 mm 

The best results in the temperature overshoot and heat flux drops are 
achieved at this amplitude. Figures 13a and 13b reflects this decrease. At the frequency of 
0.69 Hz, the overshoot decreases to its lowest magnitude. In the non-oscillating boiling 
case the maximum temperature overshoot value is 99.7 C, at this frequency, the maximum 
temperature overshoot is read as 87.6 C with 27.7% drop. The heat flux value at which 
the onset of nucleation boiling occurs, decreases up to 32.0% from 86,300 W/m^ to 58,700 
W/m^ 

The CHF value stays almost unchanged during the experiments at this amplitude. 

c. Amplitude of 0.156 mm 

Just like the amplitude of 0.078 mm, at this amplitude temperature 
overshoot and heat flux drops are achieved, but they are not as large as the amplitude of 
0.117 mm. The results of this amplitude can be seen in Figures 14a and 14b. 
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The remarkable side of this set of experiments is CHF. As can be seen from 
Table 2b, CHF value goes up from 132,000 W/m^ to 142,000 W/m^ at the same wall 
temperature as the non-oscillating case. 


Amplitude 

(mm) 

0.078 

0.117 

0,156 

Frequency 

(Hz) 

Wire 

(C) 

Bulk 

(C) 

Heat 

Flux 

(W/m^) 

Wire 

(C) 

Bulk 

Temp 

(C) 

Heat 

Flux 

(W/m^) 

Wire 

(C) 

Bulk 

Temp 

(C) 

Heat 

Flux 

(W/m^) 

No Stroke 

89.2 

51.8 

132,000 

89.2 

51.8 

132,000 

89.2 

51.8 

132,000 

0.4 

87.8 

54.0 

122,500 

87.4 

52.5 

132,000 

88.8 

53.1 

132,000 

0.56 

88.4 

53.3 

132,000 

88,0 

54.5 

132,000 

88.8 

50.8 

132,000 

0.69 

87.9 

52.3 

122,500 

87.4 

54.7 

123,000 

89.3 

51.2 

142,000 

0.8 

88.4 

53.3 

132,000 

88.0 

52.0 

132,000 

89,7 

52.2 

142,000 

1.5 

87.8 

52,6 

123,000 

87.5 

53.5 

123,000 

88.5 

53.5 

132,000 


Table 2b. Critical Heat Flux Values for Oscillating FC-72. For No Stroke Case the 
Results of Experiment 18 is Used. 


C. FURTHER DISCUSSIONS 

A total of 46 experiments have been conducted to achieve the boiling curve of FC- 
72 at various bulk temperatures, oscillation frequencies and amplitudes. The effects of the 
oscillation amplitudes and frequencies on the superheat, which is necessary to initiate 
nucleate boiling, have been observed. 

In Experiment 18 (Figure 9) without any oscillation in the bulk fluid, before the 
boiling starts the wall temperature goes up to 99.7 C and then as soon as the boiling starts 
the wire temperature drops 12.6 C. On the other hand in Experiment 39 (Figure 10) with 
the oscillation amplitude of 0.117 mm and oscillation frequency of 0.69 Hz the wall 
temperature goes up to only 87.6 C. This is the smallest temperature overshoot for 
nucleate boiling achieved in this study and it is 27.7% less than the no oscillation case 
temperature overshoot, which is Experiment 18. 
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Under all circumstances, the temperature overshoot value is reduced with 
oscillation from 8.7% (f= 0.4 Hz, A = 0.078 mm) to 27.7% (f= 0.69 Hz, A = 0.117 mm). 
The heat flux value for the onset of nucleate boiling drops in all frequencies and 
amplitudes. The drop is between 32.0% (f = 0.56 Hz, A = 0.078 mm) and 8.3% (f = 0.69 
Hz, A = 0.156 mm). 

Figure 15 shows the maximum temperature overshoot values required for boiling 
incipience in FC-72 at three different amplitudes and five different frequencies. As overall 
result at the oscillatory conditions, there is a certain decrease in required superheat. 

The amplitude of 0.117 mm is less effective at low and high frequencies, but at 
0.69 Hz the superheat reaches its lowest value. 

The reduction of superheat for the onset of boiling can be attributed to the way the 
bubbles forming at the nucleation sites depart from the hot surface. When the bubble 
completes its growth in the cavity, it starts to grow bigger on the outer surface. At this 
point in highly wetting fluids like FC-72, the bubble can not depart the surface easily and 
this occurrence results in the superheat before the boiling incipience. 

The bubble snapping process from the cavities is accelerated due to the oscillation. 
This helps to overcome the surface tension effects of highly wetting FC-72 dielectric fluid. 
As it is expected at some amplitudes like 0.117 mm, the oscillation is more effective to 
remove forming bubbles from the hot surface, but at the other two amplitudes the 
resultant reduction is not as high as 0.117 mm. This occurrence can be related to the 
active nucleation sites on the hot surface and the effective bubble size. 
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In The Following Figures from 9 to 14b These Symbols are Used to Clear the Plots: 


X Increasing Heat Flux (Wire Temp.-Saturation Temp) 
o Decreasing Heat Flux (Wire Temp.-Saturation Temp) 
+ Increasing Heat Flux (Bulk Temperature) 

Decreasing Heat Flux (Bulk Temperature) 
Kuehn-Goldstein Correlation (Natural Convection) 
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Figure 10. The Boiling Curve of FC-72 with the Lowest Temperature Overshoot. 
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Figure 12a. Top Plot in Figure 12a is the Boiling Curve of FC-72. Other Five Plots ir 
Figures 12a and 12b are the Boiling Curves with an Oscillation of Amplitude 0.078mm 
and Five Different Frequencies. 















EXP #46, AMPLITUDE=0.078mm, f=0.69Hz 
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Figure 12b. Top Plot in Figure 12a is the Boiling Curve of FC-72. Other Five Plots 
Figures 12a and 12b are the Boiling Curves with an Oscillation of Amplitude 0.078mn 
and Five Different Frequencies. 
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Figure 14a. Top Plot in Figure 14a is the Boiling Curve of FC-72. Other Five Plots ir 
Figures 14a and 14b are the Boiling Curves with an Oscillation of Amplitude 0.156mm 
and Five Different Frequencies. 







EXP #23, AMPLITUDE=0.156mni, f=0.69Hz 



T wire-T sat (C) 


Figure 14b. Top Plot in Figure 14a is the Boiling Curve of FC-72. Other Five Plots in 
Figures 14a and 14b are the Boiling Curves with an Oscillation of Amplitude 0.156mm 
and Five Different Frequencies. 
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MAX SUPERHEAT IN OSCILLATORY FC-72 



FREQUENCY (Hz) 


Figure 15. Maximum Temperature OvershootValues Required for the Onset of 
Nucleate Boiling for Three Different Amplitudes at Given Frequencies are Plotted. It can 
be Seen That at all Frequencies and Amplitudes There is a Drop in Required Superheat 
values. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The present study investigated boiling of degassed FC-72 in an oscillating fluid 
environment. An enhancement of nucleate boiling with oscillation has been confirmed by 
the results of the experiments. All of the oscillation amplitudes and frequencies changed 
the boiling curve of highly wetting dielectric fluid FC-72, so that the apparent temperature 
overshoot in non-oscillating case has decreased. Remarkably at some oscillation amplitude 
and frequencies the superheat is almost vanished. The amplitude of 0.117 mm and the 
frequency of 0.69 Hz provided the lowest wall temperature of 87.6 C with 27.7% 
decrease in the maximum temperature overshoot, which was measured to be 99.7 C for 
non-oscillating case. In this situation, the boiling incipience took place at the same wall 
temperature, at which nucleate boiling occurred. 

With the decreasing temperature overshoot, the heat flux, at which the incipience 
of boiling occurs, has dropped, too. On the other hand, critical heat flux has not been 
effected as much with the oscillation except two cases. At 0.156 mm amplitude and 0.69 
Hz and 0.8 Hz, CHF increased from 132,000 W/m^ up to 142,000 W/m^. 

The other conclusions can be drawn from this study are as follows; 

The non-oscillating boiling of degassed FC-72 gave almost identical curves at 
different runs. At different bulk temperatures, while the natural convection regime 
remained uneffected, CHF increased with decreasing bulk temperature. 

Even though the experiments, which are conducted in this study, proved to be 
reproducible within the uncertainty limits, in the nucleate boiling phase the wall 
temperature showed different behaviors as expected. 

Here it is possible to conclude that the amount of necessary temperature overshoot 
for the onset of nucleate boiling, can be decreased with oscillation in the fluid. The reason 
can be explained by the forces exerted by the oscillating fluid on the bubbles, which are 
forming on the hot wall. The oscillating fluid can detach the bubbles forming in the 
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nucleation sites as soon as they start growing on the outer surface. This increases 
efficiency of nucleation sites, which are very scarce due to the highly wetting dielectric 
fluid FC-72. It is possible that a nucleation site with increased efficiency can activate other 
nucleation sites downstream with bubble pumping effect as it has been tested by Kelleher, 
Egger, Joshi and Lloyd[Ref 8], 

The different effects of the amplitudes and frequencies tested in this study on the 
boiling curve, is related to the bubble size and growth rate, which depend on the size and 
shape of the nucleation sites. 

B. RECOMMENDATIONS 

In this study only three amplitudes and five frequencies are tested. For future 
research, a more detailed frequency and amplitude map can be created. 

With a better SEM technique, the surface of the platinum wire can be examined for 
the nucleation sites in the right range for boiling of the FC-72. 

The experimental setup can be redesign to decrease the oscillation amplitude and 
frequency uncertainties. It can be done by providing smoother piston operation for 
example by decreasing friction between the piston and the seal. 

The oscillation frequency and amplitude measurement methods can be improved 
for more accurate readings. If necessary in the future instead of oscillating the fluid, the 
same sinusoidal motion around the wire can be created by moving the platinum wire or the 
board mechanically. 

The platinum wire board with adjustable spring system can be redesigned to 
eliminate the difficulties of placing the board into the main chamber. This can be 
accomplished by providing replaceable power and measurement wires. If the wires are 
attached to the board at all times, it causes difficulties during handling. 




APPENDIX A. 


CALIBRATION OF PLATINUM WIRE 


In this study a platinum wire used as a temperature reading device. The 
objective of this calibration is to find the relationship between the temperature and the 
resistance of platinum wire. It is known that the resistance of platinum changes with 
temperature linearly and the amount of change in resistance can be calibrated to measure 
temperature 

The method, which used in this study to measure the heat flux from a surface 
and the surface temperature, has been used by Nikuyama [Ref 14] in the experiments to 
determine the boiling regime of water. He used a platinum wire to identify different 
regimes of pool boiling of water. The temperature of the wire was determined from the 
knowledge of the manner in which its electrical resistance varied with temperature. 

A. SEM ANALYSIS OF PLATINUM WIRE 

The platinum wire from the same sample, which is used during the 
experiments, is analyzed in Scanning Electron Microscope under two magnifications. The 
magnification of 925x, as it can be seen in Figure 16a, yields that the diameter of the 
platinum wire is 0.055mm. The magnification of 4930x (Figure 16b) shows the lenghtwise 
line formations, which have the width of l-2|im. Possibly these lines formed during the 
wire drawing process. At the same magnification, it is possible to see 2-3pm diameter 
bumbs on the surface. No possible nucleation sites are detected at these two 
magnification. 

B. CALIBRATION PROCEDURE IN THE CALIBRATION BATH 
During the calibration in the bath following apparatus used: 

1. Hewlett Packard 3456A Digital Voltmeter (for 4-wire ohm measurement) 

2. Rosemount Engineering Company Model 162C Serial no 985 Platinum Resistance 
Thermometer (PRT) (to measure the calibration bath temperature. According to 
Reference 15, PRT has been calibrated in comparison with a platinum resistance 
temperature standard, which is done by the National Bureau of Standards.) 
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3. Rosemount Engineering Company Model 920A Commutating Bridge (to measure PRT 
resistance for temperature measurement) 

4. Rosemount Engineering Company Model 913 A Calibration Bath (fluid in the bath is 
ethylene glycol) 

5. Rosemount Engineering Company Model 923B Power Supply (to control bath 
temperature with Neslab Endocal refrigerated circulation bath) 

The 0.05 mm platinum wire has been mounted on a plexiglass board and some 
tension applied by the spring assembly in order to keep the wire straight against thermal 
expansion. The board immersed into the Rosemount Engineering constant temperature 
liquid bath. The platinum wire electrically connected to the HP 3456A Digital Voltmeter 
in the 4-wire D measurement arrangement. The bath temperature increased from 
approximately 25 C to 85 C with the increments of 5 C. When the bath temperature is 
settled at each setting, the resistance of the platinum wire and the bath temperature 
recorded. When 85 C is reached, temperature of the bath is decreased down to 25 C with 
steps of 10 C, in order to detect possible hysteresis in measurements. 

The position of the platinum wire in the bath is very important for resistance 
readings. In the beginning of the calibration, the board placed very close to the surface of 
the bath and the resistance values on the HP 3456A display were not constant. The board, 
which holds the platinum wire, immersed into the mid sections of the bath, where the 
temperature is relatively stable and the platinum wire was next to PRT. The uncertainty of 
the resistance readings was ±0.00in. 

C. CALIBRATION PROCEDURE IN THE MAIN CHAMBER 

After the platinum wire is placed in the main chamber and the chamber is filled 
with FC-72, another calibration is conducted to check if there is any change in calibration. 
For the resistance measurements, the platinum wire connected to the HP 3456A Digital 
Voltmeter in the 4-wire SI measurement arrangement. The fluid temperature increased 
from room temperature to the saturation temperature, which is 56 C for FC-72, by the 
bulk heaters at the bottom of the main chamber. The temperature of the fluid is read by 



four thermocouples, which are mounted inside the main chamber by using HP 3497A Data 
Acquisition/ Control Unit and Unitek desktop computer with a QBASIC code. Randomly 
the resistance of the platinum wire and the FC-72 temperature recorded. 

D. RESULTS AND ERROR ANALYSIS 

The calibration data, which is conducted in the calibration bath, can be seen in 
the first two columns of Table 3 and the top plot in Figure 17. The curve fitting of 
acquired data resulted the following two linear equations; 

7” = 50.4159 ■/?-252.3321 (Al) 

7? = 0.0198 ■7’+5.0052 (A2) 

Where: T = Surface Temperature of the Platinum Wire in C. 

R = Resistance of the Platinum Wire in Q 



Table 3. Calibration Data and Linear Aggression Results. 
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The linear line in the bottom plot of Figure 17 is this curve fit line. It can be 
seen how well it matches with the calibration data points. The third column in Table 3 is 
the temperature values from Equation AI corresponding to the resistance values in the 
first column. There is error between the real temperature values and the results of linear 
aggression. The fourth column in Table 3 is the percent error between the experimental 
and curve fit temperature values. 

In-place calibration check is performed to control the changes in the calibration 
equation. Tables 4a and 4b are the results of two in-place calibrations. In Figure 18 the 
two plots are the results of these two in-place calibrations plotted with Equation Al, the 
calibration curve. 

1. The overall average of percent error in temperature calibration is 0.5327%. 
The average of nominal error between the two temperature values, namely columns 2 and 
3, is 0.2056. 

2. It is found that the correlation coefficient ‘r^’ is 0.9998 for Equation Al. 

3. The value of =0.0145 indicates that the vertical standard 

deviation of the data (from precision or random error) is only about 1.45% of the total 
vertical variation caused by the straight-line relationship between temperature and 
resistance. Namely, the ratio of the vertical standard deviation of the data about the line to 
the total vertical variation of data is 1.45%. 
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Figure 16b. SEM Picture of Sample Platinum Wire (4930x) 








TEMPERATURE (deg C) 

LINEAR LINE FITS BEST TO THE DATA 



Figure 17. 


Top Figure is the Data Acquired During the Calibration of Platinum Wire. 
The Bottom Figure Shows How the Curve Fit Matches to the 
Experimental Data. 





IN-PLACE CALIBRATION DATA 




Figure 18. The Data from Two Different In-Place Calibration Comparered 
the Curve Fit. 
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TEMPERATURE 

(degC) 

RESISTANCE 

(ohms) 

24.2450 

5.4900 

24.1350 

5.4880 

31.1400 

5.6280 

31.4800 

5.6340 

33.7500 

5.6840 

35.1400 

5.7120 

36.1600 

5.7310 

37.5200 

5.7600 

41.5100 

5.8350 

44.4800 

5.8940 

50.6700 

6.0180 

51.7700 

6.0410 

53.3200 

6.0720 

56.2200 

6.1300 

55.4900 

6.1150 

51.2400 

6.0320 


Table 4a. In-Place Calibration Data #1 


TEMPERATURE 

RESISTANCE 

(degC) 

(ohms) 

27.9040 

5.5720 

28.4540 

5.5800 

28.6160 

5.5890 

29.9420 

5.6130 

33.7830 

5.6890 

37.7440 

5.7620 

42.0970 

5.8460 

47.7220 

5.9580 

50.2260 

6.0080 

52.6040 

6.0540 

53.1780 

6.0660 

53.8960 

6.0780 

56.2070 

6.1250 

56.2240 

6.1240 

56.2300 

6.1240 

56.1840 

6.1230 

56.1950 

6.1230 


Table 4b. In-Place Calibration Data #2 


54 





APPENDIX B. SAMPLE CALCULATIONS 

A. DETERMINATION OF EXPERIMENTAL VALUES 

The following sample calculations are the results of 18th experiment, which is at 
the same time the boiling curve of the FC-72. 

1. Condenser temperature: 

The condenser temperature is the average of four thermocouple readings, which 
are embedded in the condenser surface. 

(t; + t; + r, + t,) 

cond ^ 

T’w =26.519°C 


2. Fluid bulk temperature: 

The fluid temperature is the average of four thermocouple readings, which are in 
the FC-72 medium inside the main chamber 

^ _(T,+T,+T,+T,) 

hull, - ^ 

= 53.485° C 


3. Current of platinum wire and precision resistor: 

The current of platinum wire and precision resistor is obtained from the division of 
voltage drop across the precision resistor divided by its resistance value. 





0.461 

2.02 


= = 0.2282A 



4. Platinum wire resistance value: 

The platinum wire resistance is obtained from the division of voltage drop across 
the platinum wire by platinum wire current. 


1.443 

0.2282 


= 6.32340 


5. Platinum wire surface temperature: 

To find the platinum wire surface temperature, the calibration formula in the 
APPENDIX A is used with platinum wire resistance as an input. 

= 50.4159 • - 252.3321 

=50.4159-6.3234-252.3321 
= 66.4678“ C 

6. Heat flux from platinum wire: 

The heat flux from the platinum wire is equal to the division of electrical power 
supplied to the platinum wire by the platinum wire surface area. 



„ _ 02282 • 1.443 
^ " 1.76025-10-’ 
qr" = 18,707.15W/m' 
where = tt- - ^p,»-,„ 

=^-0-1018736-5.5-10'' 
^p,,p,„ = 1.76025-10-'m' 
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7. Piston displacement: 

The piston displacement is equal to the piston cross sectional area multiplied by the 
piston stroke length. The following calculations in parts 7 and 8 are done for the stroke 
length of 1,5 cm. 


8. Amplitude of oscillation in the mair 

The amplitude of the oscillation in the n 
of piston displacement by the main chamber cro: 


n chamber is obtained from the divis: 


;r length and width of the main chamber respectively. 
0.06557 

“ 5.996-2.3377 • 


9. Frequency calculation: 

The frequency of the oscillation is the arithmetical inverse of oscillation period 
reading in seconds. The following calculations are done for lowest frequency, which 


/ = 0.4Hz 



B. DETERMINATION OF NUMERICAL VALUES 

In this part, natural convection heat transfer around the platinum wire immersed 
into FC-72 is calculated numerically by using recommended correlation by Kuehn and 
Goldstein [Ref 12] and Egger’s results [Ref 13] are used for fluid properties. 


1. Film Temperature: 


^ (53.4293 + 66.4327) 

V - 2. 

7;,„ = 59.9310° C 


2. Thermal conductivity: 


k = 


(0.6033-0.00115- 

io 


I: = 0.0534W/m°C 

3. Liquid density: 

p, =(1.740 - 0.00261 • ) • 1000 

p, = 1.5836-10’kg/m’ 

4. Kinematic viscosity: 


+273.15 
!10-®m'/s 


C, = (0.241111 + 3.70337 • 10^ ■ 7),„ ) • 4 
C, = 1.1022-10'j/kg°C 
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6. Thermal expansion coefllcient: 


7. Thermal dilTusivity: 


8. Prandit number: 


9. Rayleigh number: 


(1.740 - 0.00261 
^ = 0.001648 1/°C 

k 

a = 0.3062 10-’mVs 

Pr = - 
Pr = 9.4339 

g Op,.r,. 


10. Nusselt number: 


log 1 


2 



11. Heat Flux 


Nud = 1,4036 


D 

=19,5low/m" 
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APPENDIX C. 


UNCERTAINTY ANALYSIS 


A. 


UNCERTAINTY IN SURFACE AREA 


Dp, -SSim = 0.055/w/w 
= 4.01 \in = ]01.88/»»i 
tu P = 0002mm 

= 0.06/rt = 1.524/w/w 

= ^0.10188-5.510-‘ 
= I.76035-10-'m' 


Wire Diameter 
Wire Length 

Uncertainty in Wire Diameter 
Uncertainty in Wire Length 
Wire Surface Area 


dAp, 
S Dp, 
d Ap, 
dtp, 


= ■ Dp,,^,^ 


Then the uncertainty in surface area is; 

+{jr-Dp„ 





If 0.002^' ^ f 0 . 0 ^ 

^p,.mre \y0.0552 U.OllJ 


^ = 0.0393 or 3.93% 

Ap,W.re 

f^A.p,mre = 6.9217-10-' 
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B. UNCERTAINTY IN POWER 

2Q = 0.006% Uncertainty in Pre. Res. Vol. Drop [Ref. 11 pp 646] 

ft);, p, = 0.006% Uncertainty in Pt.Wire Voltage Drop [Ref 11 pp 646] 
ft)p = 0.020 Uncertainty in Precision Resistor 



(Vp,..] 

^ I ■R20 ) 

Then the uncertainty in power is: 



(Oq _ If 0.006%^ ^ f 0.006%^ ^ f 0.02^^ 
Q 100 J 100 J ^l 2 . 02 J 


^ = 0.099 or 0.99% 



c. 


UNCERTAINTY IN HEAT FLUX 



dA y A^J 


Then the uncertainty in heat flux is: 



^ = -J(0.0393)^ +(0.0099)^ 
^ = 0.0405 or 4.05% 

q" 


D. UNCERTAINTY IN TEMPERATURE 

The uncertainty in the thermocouple temperature measurements is the sum of the 
uncertainties due to the hardware and software, namely HP 3497A and the software used 
in this study. According to Reference 11, page 650 the temperature measurement accuracy 
with HP 3497A Data Acquisition/ Control Unit with the HP 3054 A/C Data Acquisition/ 
Control System Software is ±0.4 C in the temperature range from 0 to 300 C. Even 
though the software used in this study is different, the accuracy value from the reference 
book is good enough 
0}^.^ =0.4C 
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E. UNCERTAINTY IN WIRE SURFACE TEMPERATURE 

As it is given in Appendix A, Section C ‘Results and Error Analysis’ maximum 
uncertainty is: 

^^^ = 3,031% 

T^rf 

■=Q>lH'iC For the surface temperature of 23.27C. 

F. UNCERTAINTY IN OSCILLATION AMPLITUDE 

= 0.376/>7 

- 10">/W 

a = 5.996//J 
h = 2.3377W 
COD.Pislon = 0.005//7 

L.Piswn - 0-5W/W 

ty„ =0.01/>J 

tUj =0,01/« 
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A, Amplitude ^ DPiston LSlroke ] 

'^Amplitude * 


A, Amplitude If 0.005 ^ ^0.5^ ^ f 0.01 ^ f 0.01 

'^Amplitude ”r2-0.376J "^U.QgeJ ^13^11) 


= 0.0506 0.5,06% 

Amplitude 

G. UNCERTAINTY IN OSCILLATION FREQUENCY 

Because of the friction factor uncertainty in the frequency is higher at the lower 
frequencies. 

Tp,r,a^ = 2i5 
(»T.P.nod = 005i 

Period 

^ f _^T, Period 
f ^Period 
0 ) , 

^ = 0.02 or 2% 

/ 
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APPENDIX D. 


COMPUTER PROGRAMS 


A. MAIN DATA ACQUISITION PROGRAM 

This QBASIC program executes these operations; To set up communications with HP 
3497A and HP 6289A. To control power settings incrementally by the help of HP 59501B Power 
Supply Programmer. To keep time intervals between steps. To read thermocouple voltages and 
convert them to temperature values. To make calculations to figure out platinum wire temperature 
and heat flux. To store temperature and heat flux values in two separate data files. 


'DATA COLLECTION PROGRAM DURING THE EXPERIMENT 

■SET UP COMMUNICATIONS 

OPEN "GPIBO" FOR OUTPUT AS #1 
OPEN "GPfflO" FOR INPUT AS #2 
PRINT #1, "GPIBEOS OUT CR LF" 

PRINT#!, "ABORT" 

PRINT#!, "RESET" 

PRINT#!, "CLEAR" 

PRINT #1, "REMOTE" 

STARTIT: 

FILEIT: 

PRINT " IF A DATA FILE ALREADY EXISTS WITH THE NAME YOU SELECT, " 
PRINT " PREVIOUS DATA WILL BE OVERWRITTEN " 

INPUT " ENTER THE NAME OF THE TEMP DATA FILE TO CREATE ", FILE! $ 

IF FILE!$ = "" THEN GOTO FILEIT 
OPEN FILEIS FOR OUTPUT AS #3 

INPUT " ENTER THE NAME OF IHE FLUX DATA FILE TO CREATE ", FILE2$ 

IF FILE2$ = "" THEN GOTO FILEIT 
OPEN FILE2$ FOR OUTPUT AS #4 

'BE SURE AT THIS POINT THAT THE FLUID IS DEGASSED AND 
■CLOSE TO THE SATURATION TEMPERATURE 

CLS 

VOLTS = 20 

FOR V = ! TO !8 STEP 1 
IF V < 7 THEN 

VOLTS = VOLTS + 5 
ELSE 

VOLTS = VOLTS + 2.5 
END IF 

VALUE = 1NT(VOLTS * 10) 

VALUES = STR$(VALUE) 

IMAGES = "000" + VALUES 'CREATE LEADING ZEROS 

VOLTSS = RIGHTS(1MAGES, 3) 'SEND 3 RIGHT lUSTIFIED DIGITS 
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'SEE WHAT'S GOING OUT 


PRINT "COMMAND SENT WAS VOLTSS 
PRINT "" 

PRINT #1, "OUTPUT 6;2" + VOLTSS 
TIMEPAST = 0 
STARTTIME = TIMER 
WHILE TIMEPAST < 60 

TIMEPAST = TIMER - STARTTIME 
WEND 

FINISHS = TIMES 
FORP= 1 TO 10 

' TIME DELAY BETWEEN DATA POINTS 
TIMEPAST = 0 
STARTTIME = TIMER 
WHILE TIMEPAST < 3 'DELAY TIME IN SECONDS 

TIMEPAST = TIMER - STARTTIME 
WEND 

FINISHS = TIMES 

' READ VOLTAGES ON PRECISION RESISTOR AND PLATINUM WIRE 

PRINT # 1, "OUTPUT 9 ; AR AF20 AL21" 'READ CHANNELS FROM 20 TO 21 

FOR I = 0 TO 1 

PRINT #1, "OUTPUT 9; AS" 

PRINT #1, "ENTER 9" 

INPUT #2, DATS 
VOLT(l) = VAL(DATS) 

NEXT I 

' CALCULATIONS 

VOLT20 = VOLT(O) 

VOLT21 = VOLT(l) 

AMP20 = VOLT20 / 2.02 
AMP21 = AMP20 
RPLAT = VOLT2I / AMP21 
TPLAT = 50.4159 RPLAT - 252.3321 
QFLUX = VOLT21 ■* AMP21 / .000017602529# 

TOTALV = VOLT(0) + VOLT(l) 

' COLLECT DATA FOR BULK, CONDENSER AND AMBIENT TEMPERATURES 
PRINT #1, "OUTPUT 9; AR AFOO AL08" 'CHANNELS 00 THRU 08 

FOR I = 0 TO 8 

PRINT # 1, "OUTPUT 9;AS" 'ANALOG STEP 

PRINT #1, "ENTER 9" 

INPUT #2, DATS 

VOL(I) = VALIDATS) 'CONVERT STRING TO NUMBER 

E = VOL(I) 

C(l) = .10086091# ■►E''0 + 25727.94369#'"E^ 1 -767345.8295# ■* E''2 + 78025595.81# • E3 
- 9247486589# » E ^ 4 + 697688000000# * E'' 5 - 26619200000000# » E '' 6 + 39407800000000# E '' 


'ANALOG STEP 

'CONVERT STRING TO NUMBER 

PRECISION RESISTOR 
'PLATINUM WIRE 
'PRECISION RESIS.= 2.02 ohms 

'AREA=1.7602529E-5 sq.meter 


'"2" SETS HIGH RANGE (10 VOLTS) 

'DELAY TIME IN SECONDS 



NEXTl 

AMBI = C(8) 

COND = (C(l) + C(2) + C(3) + C(0)) / 4 
BULK = (C(5) + C(6) + C(7) + C(4)) / 4 
PRINT ■' BULK TEMP. = "; BULK 

PRINT " CONDENSER TEMP. = COND 
PRINT 

DELTAT = TPLAT - BULK 
' PRINT DATA TO SCREEN 

PRINT 'TIME: FINISHS; " VOLTS= VOLTS; 

PRINT " V= V; " P= P 

PRINT "TPLAT = ", TPLAT; " DELTA T= "; DELTAT; " TOTAL VOLT = "; TOTALV; " "; 
VOLT21 

PRINT "QFLUX= "; QFLUX 
PRINT "" 

' SAVE DATA TO FILE 

PRINT #3, USING "####.###"; BULK; COND; AMBI; DELTAT; TOTALV; VOLT21 
PRINT #4, USING "######.###"; QFLUX 
NEXTP 

NEXT V 

FOR V = 17 TO I STEP -I 
IF V > 5 THEN 

VOLTS = VOLTS - 2.5 
ELSE 

VOLTS = VOLTS - 5 
END IF 

VALUE = INT(VOLTS ♦ 10) 

VALUES = STR$(VALUE) 

IMAGES = "000" + VALUES 
VOLTSS = RIGHTS(IMAGES, 3) 

PRINT "COMMAND SENT WAS "; VOLTSS 
PRINT "" 

PRINT #1, "OUTPUT 6;2" + VOLTSS 
TIMEPAST = 0 
STARTIIME = TIMER 
WHILE TIMEPAST < 60 

TIMEPAST = TIMER - STARTTIME 
WEND 

FINISHS » TIMES 
FORP= 1 TO 10 

' TIME DELAY BETWEEN DATA POINTS 
TIMEPAST = 0 
STARTTIME = TIMER 
WHILE TIMEPAST < 
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■CREATE LEADING ZEROS 
■SEND 3 RIGHT lUSTIFIED DIGITS 
'SEE WHAT'S GOING OUT 

'"2" SETS HIGH RANGE (10 VOLTS) 


'DELAY TIME IN SECONDS 


'TIME DELAY IN SECONDS 




TIMEPAST = TIMER - STARTTIME 
WEND 

FINISHS = TIMES 

' READ VOLTAGES ON PRECISION RESISTOR AND PLATINUM WIRE 

PRINT #1, "OUTPUT 9 ; AR AF20 AL21" 'READ CHANNELS FROM 20 TO 21 

FOR I = 0 TO 1 

PRINT #1, "OUTPUT 9; AS" 

PRINT #1, "ENTER 9" 

INPUT #2, DATS 
VOLT(I) = VALPATS) 

NEXT I 

’ CALCULATIONS 

VOLT20 = VOLT(O) 

VOLT21 = VOLT(l) 

AMP20 = VOLT20 / 2.02 
AMP21 = AMP20 
RPL AT = VOLT21 / AMP21 
TPLAT = 50.4159 * RPLAT - 252.3321 
QFLUX = VOLT21 * AMP21 / .000017602529# 

TOTALV = VOLT(O) + VOLT(l) 

' COLLECT DATA FOR BULK, CONDENSER AND AMBIENT TEMPERATURES 
PRINT #1, "OUTPUT 9; AR AFOO AL08" 'CHANNELS 00 THRU 08 

FOR I = 0 TO 8 

PRINT # 1, "OUTPUT 9;AS" 'ANALOG STEP 

PRINT #1, "ENTER 9" 

INPUT #2, DATS 

VOL(I) = VALPATS) 'CONVERT STRING TO NUMBER 

E = VOL(I) 

C(I) = . 10086091# * E-'0 + 25727.94369# » E'' 1 - 767345.8295# » E'' 2 + 78025595.81# * E ^ 3 
- 9247486589# * E 4 + 697688000000# * E '' 5 - 26619200000000# E '' 6 + 39407800000000# E '' 
7 

NEXT I 
AMBI = C(8) 

COND = (C(l) + C(2) + C(3) + C(0)) / 4 
BULK = (C(5) + C(6) + C(7) + C(4)) / 4 
PRINT " BULK TEMP. = BULK 
PRINT " CONDENSER TEMP. = COND 
PRINT "" 

DELTAT = TPLAT - BULK 

' PRINT DATA TO SCREEN 

PRINT "TIME: FINISHS; " VOLTS= VOLTS; 

PRINT" V= ";V;" P= ”; P 

PRINT "TPLAT = TPLAT; " DELTA T= ”; DELTAT; " TOTAL VOLT = TOTALV; " "; 
VOLT21 

PRINT "QFLUX= "; QFLUX 


'ANALOG STEP 

'CONVERT STRING TO NUMBER 

'PRECISION RESISTOR 
'PLATINUM WIRE 
'PRECISION RESIS.= 2,02 ohms 

'AREA=1.7602529E-5 sq.meter 



PRINT"" 


' SAVE DATA TO FILE 

PRINT #3. USING "####.###"; BULK; COND; AMBI; DELTAT; TOTALV; VOLT21 
PRINT #4, USING "######.###"; QFLUX 

NEXTP 

NEXT V 

CLOSE #3 
CLOSE #4 

QUES: 

INPUT " DO YOU WISH TO TRY ANOTHER RUN ( Y/N ) ? ", AGAINS 
IF AGAINS = "Y" THEN GOTO STARTIT 
IF AGAINS = "y" THEN GOTO STARTIT 
IF AGAINS = "N" THEN GOTO ENDIT 
ff AGAINS = "n" THEN GOTO ENDIT 
GOTO QUES 
ENDIT: 

PRINT # I, "CLEAR " 'CLEAR ALL INSTRUMENTS ON BUS 

PRINT #1, "LOCAL " 'PLACE ALL INSTRUMENTS IN LOCAL MODE 

END 


B. DATA EVALUATION AND PLOTTING CODES 

The data stored in two files are plotted by these three MATLAB codes (deney 18,m, 
average.m, kuehn.m) with Kuehn and Goldstein correlation results. Egger’s FORTRAN code [Ref. 
13, Appendix D] has been used as a model for correlation calculation. 

a. Matlab Code ‘DENEYI8.M’ 

% THIS PROGRAM IS USED IN THE THESIS 'BOILING OF HIGHLY WETTING LIQUIDS IN 
% OSCILLATORY FLOW' TO PROCESS AND PLOT EXPERIMENTAL DATA WITH 
% FUNCTIONS average(temp,(lux) AND kuehn(tbulk,twire) 

clear 

load temp 18 
load flux 18 
temp=tempI8; 
flu.x=fluxI8; 


|f,tbulk,twire]=average(temp,flux), 

q=kuehn(tbulk,twire); 



dt=twi re-56; 


semilogy(dt,f,dl(l:18),f(l:18),'x',dt(19;35),fl:i9;35),'o',dt(l:8,:),q(l;8,;),tbidk(l:18),f(l:18),'+',tbulk(19:35) 

,f(19:35),'.');grid 

ylabelCHEAT FLUX (W/sq.m)') 

xJabelCT wire - T sat (C)’) 

tiUeCEXP #18, AMPLITUDE=Oinm') 

axis([0 70 10000 1000000]) 

print deney 18 -deps 

b. Matlab Code ‘AVERAGE.M’ 

% THIS FUNCTION IS USED IN THE THESIS ’BOE-ING OF HIGHLY WETTING LIQUIDS IN 
% OSCILLATORY FLOW’ TO AVERAGE EXPERIMENTAL DATA TO BE USED IN 
% MAIN MATLAB CODE deney_.m WITH FUNCTION kuehn(tbulk,twire) 


funclion[f,tbulk,twirel=average(temp,flux) 

fori=l :35 
ttotal=0; 
ftotal=0; 
bulktotal=0; 
wi retotal=0; 

forj=l:10 

b=(((i-l)*10)+j); 
bulktotal=temp(b, 1 )+bulktotal; 
wiretotal=temp(b,4)+temp(b,l)+wiretotal; 
ftotal =flux(b)+ftotal; 
end 

tbulk(i)=bulktotal/10; 
twire(i)=wiretotal/10; 
f(i)=ftotal/10; 

fHT; 

tbulk=tbulk'; 
twire=twire'; 


% heat flux 
% bulk temp. 

% wire surface temp. 


c. Matlab Code ‘KUEHN.M’ 

% THIS FUNCTION IS USED IN THE THESIS 'BOILING OF HIGHLY WETTING LIQUIDS IN 
% OSCILLATORY FLOW TO CALCULATE NATURAL CONVECTION HEAT TRANSFER 
% AROUND PLATINUM WIRE WITH THE CORRELATION RECOMMENDED BY KUEHN 
% AND GOLDSTEIN. 

% MAIN MATLAB CODE deney_.m USES THIS FUNCTION WITH average(temp,flux) 


function[q]=kuehn(tbulk,twire) 



% diameter of the wire (m) 


d=50e-6; 

% gravitational acceleration (m / sec''2) 
g=9.81; 

% bulk temperature (C) 
tbulk; 

% film temperature (C) 
tfilm=(tbulk+twire)/2; 

% thermal conductivity (watts/ m , K) 
k=(.6033-.00115*tfilm)/10; 

% thermal expansion coefficient (1 / K) 
beta= 0026l./(1.74 - .00261»tftlm); 

% kinematic viscosity of the fluid (sq.m / sec) 
nu= 1.203952e-8*exp( 1058.4109./(tfilm+273.15)); 
% specific heat of the fluid (J / kg . K) 
cp=(.2411 l+3.70374e-4*tfilm)»4186; 

% liquid density of FC-72 (kg / cu.m) 
ro=( 1.740-.00261 •tfilm)* 1000; 

% thermal diffusivity of the liquid (sq.m / sec) 
alfa=k./ro./cp; 

% prandlt number of the fluid 
Pr= nu ./ alfa; 

% rayleigh number of the fluid 
Ra=g*beta.*(twire-tbulk).*d''3./(nu.*alfa); 

% calculate nusselt number 


terl=(.559./Pr).''(3/5); 






lerm=( (.518*Ra.''.25.*(l+terl).''(-5/12)) .''IS + (.l*Ra.''(l/3)) .''15 

Nu=2./(log(l+(2./term))); 

% calculate heat flux (watts / sq.m) 
q=k.*Nu.*(twire-tbulk)/d; 
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